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Abstract 

Here I give a more detailed account of the part of the conference reportQ that was devoted 
to reinterpreting the Einstein 'unified models of gravity and electromagnetism' (1923) as the 
unified theory of dark energy (cosmological constant) and dark matter (neutral massive vector 
particle having only gravitational interactions). After summarizing Einstein's work and related 
earlier work of Weyl and Eddington, I present an approach to finding spherically symmetric 
solutions of the simplest variant of the Einstein models that was earlier mentioned in Weyl's 
work as an example of his generalization of general relativity. The spherically symmetric static 
solutions and homogeneous cosmological models are considered in some detail. As the theory 
is not integrable we study approximate solutions. In the static case, we show that there may 
exist two horizons and derive solutions near horizons. In cosmology, we propose to study the 
corresponding expansions of possible solutions near the origin and derive these expansions in 
a simplified model neglecting anisotropy. The structure of the solutions seems to hint at a 
possibility of an inflation mechanism that does not require adding scalar fields. 

1 Introduction 

In this report I give a new interpretation of the 'unified theory of gravity and electromagnetism' 
proposed by A.Einstein in 1923 in pQ and briefly summarized in [2j. Einstein gave no details of his 
derivations, presented no exact or approximate solutions, and did not explain why he completely 
abandoned his theory (I failed to find any reference to his papers [T] - [2] in his later work). 
Apparently these papers were soon forgotten by the scientific community and I could not find any 
reference to these papers in the second half of the 20-th century except for interesting remarks by 
Schrodinger [3] and a critical discussion by Pauli in addenda to the English translation of his famous 
book [4|. For these reasons, I first give a brief historical introduction summarizing Einstein's ideas 
and results as well as earlier related work of Weyl and Eddington. 

Immediately after the general relativity was formulated in its final form (1915 -1916) some 
attempts to modify it started. Einstein himself added the cosmological constant term A to save 
(unsuccessfully) his static cosmology. After Friedmann's work (1922-1924) this modification was 
becoming more and more dubious. Weyl, after 1918, developed a much more serious modification 
aimed at unifying gravity and electromagnetism (most clearly summarized in [5]). Starting from 
Levy-Civita's ideas on a general (non-Riemannian) connection (1917) he developed the theory of 
a special space in which the connection depends both on metric tensor and on a vector field which 
he tried to identify with the electro-magnetic potential. To get a consistent theory he introduced 
a general idea of gauge invariance which survived although the theory itself failed as he admitted 
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later. In paper [6 J Einstein discussed Weyl's theory and expressed (like Pauli in [1]) the opinion 
that the theory is mathematically very interesting but probably not physical, at least, in its original 
formulation. 

In 1919 Eddington proposed a more radical modification of general relativity [7], [8]. His idea 
was to start with the pure affine formulation of the gravitation, i.e. using first the general symmetric 
affine connection and only at some later stage introducing a metric tensor. Indeed, the curvature 
tensor can be defined without metric (here we use Einstein's notation [1] but denote differentiations 
by commas): 

i _ _yi I yi -pn , -pi _ -pi j-rn fi\ 
' klm 1 kl,m ~r~ 1 nl 1 km "i - 1 km,l 1 nm 1 - kl • \ ± J 

Then the Ricci-like (but non-symmetric) curvature tensor can be defined by contracting the indices 
i, m (or, i, I): 

r kl — — 1 kl,m + 1 nl L km + 1 km,l ~ 1 nm L kl \ z ) 

(let us stress once more that T™ z = but r k i ^ r^). Using only these tensors and the anti- 
symmetric tensor density one can build up a rather rich geometric structure. In particular, Ed- 
dington discussed different sorts of tensor densities [8] . A notable scalar density is 



C = J-det(r kl ) = y/^F (3) 



which resembles the fundamental scalar density of the Riemannian geometry, \/ —det(gki) = y/—g- 
For this and some other reasons Eddington suggested to identify the symmetric part of r^i with 
the metric tensor. The anti-symmetric part, 

<f>kl = 2 (XTm,l ~ T lm,k) i <\>kl,m + 4>lm,k + 4>mk,l = , (4) 

strongly resembles the electro-magnetic field tensor and it seems natural to identify it with this 
tensor. Eddington tried to write consistent equations of the generalized theory but this problem 
was solved only by Einstein. 

The starting point of Einstein in his first paper of the series [1] was to write the action principle 
and to suppose ([3]) to be the Lagrangian density depending on 40 connection functions T^. Varying 
the action w.r.t. these functions he derived 40 equations that allowed him to find the general 
expression for T 1 ^ (the derivation is similar to that of the standard general relativity) : 

T m = l [g m» (v( + ^ _ SM n) _ SM e + 1 {6 m fj + ^ - fc)] (5) 

Here Ski is a symmetric tensor (s mn is the inverse matrix to Ski), which Einstein interpreted as the 
metric tensor (then the first term is the Christoffel symbol for this metric), and i n is a vector which 
he tried to connect with the electro-magnetic field. This identification apparently follows from the 
equations 

rkl = Rkl + ^[(ik,l ~ H,k) + hk] , (6) 

D 

<Pm = - n,k) , (7) 

which can be obtained by inserting the expression © into ©, (j3|); R k i is the standard Ricci 
curvature tensor for the metric sm ■ Einstein's interpretation of cpki a s the Maxwell field is not 
so natural because of the term ij.ii m the r.h.s. of Eq.([H]) which in fact makes this interpretation 
impossible. First, this term is not gauge invariant (but the gauge invariance was not yet discovered, 
the first clear formulation of the gauge principle was given by V.Fock in 1926). For Einstein, the 
main problem was that the electro-magnetic field in this theory could not exist without charges (i.e. 
there is no free field). To solve this problem he suggested to make this term 'infinitesimally small' 
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by choosing the corresponding dimensional constant (above, we omit all dimensional constants 
that can easily be restored). But we, today, cannot be satisfied with this solution because this 
term violates gauge invariance and makes the photon effectively massive (while it is known that 
there exist no continuous transition from the massless to massive photon theory ]q We return to 
discussing these facts, on which our interpretation of the Einstein theory is based, after considering 
the final proposals of Einstein. 

In his first paper {'Zur allgemeinen Relativitatstheorie'), he considered two limiting cases. He 
showed that, when the i n -terms in the connection vanish, the theory is equivalent to the standard 
general relativity with the cosmological term that emerges naturally and cannot be removed. In 
the flat space limit he demonstrated that weak fields (fiki (linear approximation) satisfy the free 
Maxwell equations provided that the i m z n -terms can be neglected. In the second paper of the series 
[lj he gave the following expression for the effective Lagrangian density: 



t = -2^-det(r mn ) +R- X -s mn i m i n . (8) 

This should be varied w.r.t. §ki and fki, which are the tensor densities denned with the aid 
of the scalar density \J — det(ski) and corresponding to the tensors in the decomposition, 

Tki = su + 4>ki ; (9) 

R is the scalar curvature density for the metric sjy. The Lagrangian (|8|) contains a very complex 
term \J — detr mn which is more general than the so called Born-Infeld Lagrangian proposed ten 
years later [10] (the first attempts to construct nonlinear electro-dynamics were undertaken in 
|llj). Apparently, Einstein did not try to find any particular solution of this theory and, instead, in 
the beautiful third paper ""Zur affinen Feldtheorie 1 (pQ) he proposed a significantly simpler effective 
Lagrangian that is the main subject of this paper. As he mentioned in the first two papers the 
actual form of the Lagrangian is unimportant for getting the connection the only important 
thing is on which variables it depends. 

The main idea of the third paper is to take for the Lagrangian C an arbitrary function of Ski and 
4>ki H Then he introduces the Legendre transformation and the transformed (effective) Lagrangian 
density £*: 

= 9 H , Tr- = f kl ; a«=«xn, <t>hl = ^ (10) 



Introducing the Riemann metric tensor g k i and the i fc -vector, 

g kl V=g = g kl , gu g lm = C ; i k = dj kl , (11) 

he claims (without proof) that Eq.([S]) is valid with s^i replaced by g k i and thus the affine geometry 
is the same for any £(ski, <t>kl)- Finally, he uses the freedom in choosing C*(g kl ,f kl ) and proposes 
the following effective Lagrangian density: 

C* = 2aV^g-l(3fkif kl , (12) 

where a and (3 are some constants not defined by the theory. This Lagrangian incorporates main 
properties of the theory discussed in two previous papers but is easier to deal with. 



2 The present best experimental upper bound on the photon mass is m 1 < 10~ 51 g [S]- Theoretical wisdom says 
that it must be zero. 

3 In his previous work Einstein implied that C depends on r^i . At this point he quoted an unpublished work of 
'Droste (from Leiden)' who 'two years ago expressed similar views'. The meaning of this rather cryptic remark is 
clarified in the second brief paper of ref. [2], where he confirms that (Johannes) Droste proposed to use a similar 
effective Lagrangian and, possibly a similar model (maybe, without cosmological constant). 
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To further clarify the relation of the new theory to general relativity Einstein rewrites the 
Lagrangian so that the equations of motion can be obtained by varying it in the metric and the 
vector field tensors, g k i and f k i ■ Neglecting dimensions (for example, taking h = c = k = 1) and 
changing Einstein's notation we write it as follows: 

t = V=9 [R ~ 2A - F kl F kl - n 2 A k A% F kl = A k>l - A k>l . (13) 

Now it is absolutely evident that the vector field A k is not the Maxwell field Obviously, A k is a 
neutral massive vector field with coupling to gravity only. We will call it vecton, that is an old 
fashioned but proper term for this 'geometric' particle. This particle has not been directly observed 
but it can be considered as one of the possible candidates for dark matter. In view of the fact that 
the affine theory also predicted the cosmological constant term which is one of the best candidates 
for explaining dark energy, Einstein's theory may be considered as the first unified model of 
dark energy and dark matter. 

Before we turn to further study of this model let us finish our presentation of its history. If 
you compare the Einstein model with the concrete models proposed in Weyl's book [S], you will 
find that Lagrangian similar to Eq. (ll3p is one of Weyl's examples. Einstein's and Weyl motivations 
and approaches were quite different, and the Weyl connection does not coincide with Eq.([2]) (see 
Addendum). Weyl's approach was mostly geometrical and he wrote the Lagrangian as a simplest 
illustration of possible physical applications, responding to criticism by Einstein, Pauli and other 
physicists. Einstein was most interested in physics and, especially, in cosmology. Weyl criticized 
Einstein for his departure from geometric foundations of physics, in particular, for his derivation 
of geometry from the variational (action) principle which, probably, was his main achievement in 
the third paper. Note also that Weyl included the cosmological term only to avoid contradiction to 
Einstein cosmology of that time ('before Friedmann') while in the original Einstein model ([3]), ([7]) it 
was unavoidable. I think that, conceptually, the model (I13p is a step backwards, in comparison with 
the original theory, ©, ©. There were, probably, two reasons for this step. First, Einstein's deep 
belief in simplicity of fundamental laws ('...aber boshaft ist Er nicht 1 ). Second, his disappointment 
in static cosmology after accepting Friedmann's results, [12j . Anyway, in his last papers on affine 
theory [2] he set the cosmological term to zero what is impossible in the original theory and quite 
unnatural in the framework of the affine approach. 

Above, we also mentioned work and ideas of Eddington. The intensive exchange of ideas 
between Einstein, Weyl and Eddington resulted in interrelations in their work (published in 1918- 
1923) that are difficult (and, possibly, unnecessary) to disentangle. As the constructive ideas 
of the affine theory were mostly created by Weyl, Eddington, and Einstein, the resulting model 
should probably be called Einstein- Weyl-Eddington unified model of dark energy and 
dark matter. However, as far as I am here discussing the concrete Lagrangian (|13p . I call it 
Einstein- Weyl model. 

Before turning to new results let us briefly summarize the results and thoughts of Weyl, Edding- 
ton, and Einstein. 1. Weyl had a very clear and original geometric ideas, but: a) his physics was 
rightly criticized by Einstein, Pauli, and other physicists, b) he considered the theory as a unified 
theory of gravity and electromagnetism but his vector field was also not electro-magnetic, c) his 
discussion of dynamics was incomplete and he himself regarded it as preliminary. Nevertheless, 
it is possible that not all the potential of the Weyl ideas is understood and used. 2. Edding- 
ton proposed to use, instead of the Weyl's non-Riemannian 'metrical spaces', the most general 

4 Einstein tried to identify Ak with a 'cosmic current' (this explains his notation ik). A similar identification 
reappeared much later in quantum field theory (in the vector dominance model) under the name 'field-current 
identity'. However, it is meaningless in the classical Maxwell theory. 

5 In 1917 de Sitter discovered a non-static solution of the empty - space Einstein equations with the cosmological 
constant. In 1923 Weyl and Eddington found the effect of recession of test particles in the de Sitter universe. Thus 
A was becoming useless and Einstein finally dismissed it in 1931, after Hubble's discovery (1929). 
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spaces with symmetric affine connection (without torsion). He discussed possible invariants that 
can be used in physics, in particular, the square root of the determinant of R^i EH He proposed to 
consider the symmetric part of the curvature matrix as the metric in the general space and the 
anti-symmetric one as the electro-magnetic field tensor. In later works he discussed a possibility 
to use this as a Lagrangian (long before the proposal of Born and Infeld). However, he did not 
find a consistent approach to dynamics. 3. Einstein started with formulating dynamics by use 
of the Hamilton principle similar to one proposed by Palatini in general relativity. The new (and 
crucial) idea was not to introduce any metric at the beginning and not to fix any special form of the 
affine connection (apart of the symmetry condition). He soon realized (in the second apper), that 
he does not need to use a concrete form of the Lagrangian that can be just any function (tensor 
density) of Ski and ^-matrices (see (f MTTj) ), For any such Lagrangian he proved that the affine 
connection allows one to introduce a symmetric metric and found the expression for connection. 
Both Einstein's and Weyl's expressions are special cases of the general formula for the symmetric 
connection (see Appendix). 

The most important thing is the following: supposing that the equations of motion follow from 
an action principle with the general Lagrangian fixes the geometry (connection) and, eventually, 
allows one to fix some metric compatible with this non-Riemannian connection. Another important 
thing is that the action can be (and should be) written without metric. Using the Legendre 
transformation Einstein bypassed difficulties that were met on this way and wrote more tractable 
effective Lagrangian, but some conceptually beautiful and important features of his new theory 
were thus hidden (or even lost). 

Apparently, Einstein was disappointed in the cosmological constant and also gradually realized 
that his interpretation of the anti-symmetric field as the electro-magnetic field was not quite satis- 
factory. Anyway, he completely abandoned this model and left no detailed account of his work. He 
did not mention any static or cosmological solutions even in the simplified version of the theory, 
(I13p . In this paper we try to fill this gap and establish grounds for comparing this model to the 
present day cosmology. 

2 Spherical reduction - static and cosmological solutions 
2.1 Vecton-dilaton gravity 

At first sight, the theory (|13p is very close to the well-understood Einstein-Maxwell theory which 
can be obtained when = 0. However, we will show that the two theories are qualitatively different 
and it is hardly possible to construct a reasonable perturbation theory in the parameter fj, 2 . We 
start our qualitative analysis without assuming that this parameter is small. The natural object for 
this analysis is the spherically reduced theory. When /i = 0, the theory automatically reduces to 
rather simple one-dimensional equations that can be explicitly solved. The solution is the Reissner 
- Nordstrom black hole (when the electric charge vanishes it reduces to the Schwarzschild black 
hole). In general, when gravity couples to other (not electro-magnetic) fields the spherically reduced 
theory is described by two-dimensional differential equations which are not integrable except very 
special cases (for many examples and references see, e.g., [T3]-|17j.). 

Following the approach to dimensional reduction and to resulting 1+1 dimensional dilaton 
gravity (DG) developed in papers [ISJ-pJJ] it is not difficult to derive these equations. The general 
spherically symmetric metric is = 0, 1; x° = t,x 1 = r): 



6 In addendum to [8] , where he gave a clear and detailed account of the Einstein work [T] , he also discussed another 
scalar density that was proposed by R. Weitzenbock 




(14) 
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Supposing that all other functions also depend on t, r, inserting the metric (|14p into the action 
with the Lagrangian (|13p . and integrating out the angle variables 9, (j) one can derive the following 
effective two-dimensional Lagrangiar0: 

£( 2 ) = ^[<pRW + 2 - 2Ap + {dipfjl^ -<pFijF ij - (pfjpAtA*], (15) 

where is the two-dimensional Ricci curvature depending on the gu . It is convenient to remove 
the fourth term by the Weyl rescaling of the metric, gij = (p~zgjj . Below we use the transformed 
Lagrangian, 

C<§ = ^-g[^RW + 2v- 1/2 -2kv ll2 -^ l2 F 2 -w 2 A 2 ]. (16) 

It is easy to derive the equations of motion which in a generic metric are equivalent to the 
Einstein equations for the spherically symmetric solutions of the four-dimensional theory (|13|) . By 
varying w.r.t. the diagonal metric functions ga we fist write the energy and momentum constraints. 
In the light cone (LC) metric, ds 2 = —Af(u,v)dudv, these constraints are simple: 

fdiidiv/W + tp^A 2 = 0, i = u,v. (17) 

The constraints (I17p should be derived using the general metric g^ . The other equations of motion 
may be obtained directly in the LC-metric: 

8 u d v if + f(2y- 1 ' 2 - 2Ai P 1 / 2 - l -^l 2 f~ 2 F 2 v ) = 0, F uv = A U)V - A v>u , (18) 

8 3 (^ /2 r l F lJ ) = tp^Aj i,j = u,v. (19) 

From the last equation immediately follows that d v (<pA u ) + d u ((pA v ) = . In the original four- 
dimensional theory this is the d^^—gA^) = condition eliminating spin 0. Weyl, Eddington and 
Einstein called it the Lorentz condition although we know that its origin and meaning are quite 
different from the gauge fixing condition in the Maxwell theory first introduced by L.Lorenz and 
later popularized by H. A. Lorentz. 

This dilaton gravity coupled to massive vector field (I suggest to call it vecton-dilaton gravity, 
VDG) is more complex than the well studied models of dilaton gravity coupled to scalar fields and 
thus it requires a separate study. The natural first question is: are there exact analytical solutions 
like Schwarzschild or Reissner-Nordstrom black holes? If the vector field is constant, we return 
to exactly soluble DG having explicit solutions with horizons. Otherwise, when the vector field is 
nontrivial, the answer is more difficult to find but it is worked out in some detail below. The second 
question is: what are the simplest cosmological solutions in this theory? Thus, the first thing to 
do is to further reduce the theory to static or cosmological configurations. Consider first the static 
reduction. 

2.2 Static states and horizons 

The simplest way to derive the corresponding equations is to suppose that all the functions in the 
equations depend on r = u + v . But this is not the most general dimensional reduction of the two- 
dimensional theory. There exist more general ones that allow us to simultaneously treat black holes, 
cosmologies and some waves. These generalized reductions were proposed in papers [20], [21], [T7] 
devoted to dilaton gravity coupled to scalar fields and Abelian gauge fields; here we only discuss in 
some detail the static and cosmological reductions. In both cases it can be seen that the perturbed 

7 Very similar effective Lagrangians can be obtained from the higher-dimensional analogs of the Lagrangian (|13[) . 
On the (l+l)-dimensional level it is not difficult to include other sorts of matter that appear in reductions of higher- 
dimensional supergravity theories (for references see, e.g., [T5] - |17|). I hope to discuss some of these generalizations 
in future publications. 
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theory (with a nonvanishing mass term) is qualitatively different from the non-perturbed one. 
Indeed, the non-perturbed theory is just dilaton gravity coupled to electromagnetism. This model 
is equivalent to pure dilaton gravity, which is a topological theory. In particular, it automatically 
reduces to one-dimensional static or cosmological models that can be analytically solved. Static 
states are the Reissner-Nordstrom black holes perturbed by the cosmological constant and having 
two horizons, while the space between horizons may be considered as an unrealistic cosmology. 
This object is known from 1916 times; certainly it was familiar to Einstein in 1923 but he did not 
discuss the static configuration and apparently did not consider black holes or horizons as having 
any relation to physics. 

Let us now write the static equations corresponding to the naive reduction to one spatial 
dimension. To obtain them one can reduce either the equations or the Lagrangian. Following [19J, 
|22j . we write the equations of motion in a somewhat unusual form. Let us define two additional 
functions, x an d B, by the equations (the prime denotes differentiations w.r.t. r) 

p'{r) = x(r), A'(r) = f(r)<p-V 2 (r)B(r) , (20) 

where, as follows from Eq. (|19|) . A v (r) = —A u (r) = —A(r). Then the other equations are 

X ' = -fU, B' = ~w 2 A, f = (f/ x )[-fU + w 2 A 2 ], (21) 
where we defined the potential 

U = 2(^,-1/2 - _ ^-3/2^2) ) (22) 

These equations are not integrable and cannot be solved analytically. To get numerical 
solutions we first have to study the analytic and asymptotic properties of their solutions. 

Here we only consider solutions near possible horizons that are defined as zeroes of the 
metric, / — > for finite values of <p — > (po . It is not difficult to understand that we also should 
require that A is finite near the horizon. To study the behaviour of the solutions for small values of 
(p = ip — cpo it is most convenient to consider the solutions as functions of cp. Further analysis 
shows that the solutions can be expanded in power series of <p and that the functions F = f jx and 
A = A/x should be finite. Thus we have: 

F'{p) = pF(p) » 2 A 2 (<p) , X '(p) = ~H<P) U (p) . ( 23 ) 

B'(p) = -\p^A{p) , A'(p) x (p) = Hp) lP~ 3/2 B&) + U(<p)A{<p)] , (24) 

where now the prime denotes differentiation in the new variable tp. It is not very difficult to show 
that <po, Aq, Bq, Fq can be taken arbitrary up to one relation that should be satisfied due to the 
second equation (f24l) : 

A U + p 3/2 B = 0, U = U(<po,B ). (25) 

This equation can be solved w.r.t any parameter. It is interesting to see that it has two solutions 
for tpo which means that there may exist two horizons^] as distinct from the Schwarzschild black 
hole. Note that the solutions with different Fq are equivalent because the equations are invariant 
under the scale transformation F => CF, x Cx- 

Now, following the method of [22j . one can find several terms in the expansion of the solution. 
Unfortunately, it is not clear how to construct the complete expansion and therefore our derivations 
do not allow us to study global properties of the solutions. They say nothing about asymptotic 

8 This is similar to the charged Reissner-Nordstrom black hole or to black holes in higher - dimensional theories, 
|23| . although in the present model there is no conserved electric charge. 
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properties and singularities which should be the subject of separate investigations |£| When the 
qualitative properties of the black hole type solutions will be understood, the static solutions 
and their formation can be studied by numerical simulations. As far as I know, the coupling of 
massive neutral vector particles to gravity did not attract much attention (see, however, numerical 
simulations of the critical collapse of a massive vector field in [23] ) . 

2.3 Cosmology 

2.3.1 General formulation 

Let us turn to cosmological reductions. The simplest cosmology can be obtained by the same naive 
reductions as was used for static states. However, this is not the most general dimensional reduction 
giving all possible spherically symmetric cosmological solutions (similarly, the above naive reduction 
does not give all possible static spherically symmetric solutions). A more general procedure is 
described in |21j. Following this procedure we return to the two-dimensional Lagrangian (I15D but 
add to it a scalar field ip that should represent 'matter' in low dimensions. To simplify comparison 
to the standard cosmological solutions let us now use the (t, r)-coordinates and write the general 
spherically symmetric metric as: 

ds\ = e 2a dr 2 + e 2l3 dn 2 (9, (f>) - e 2 ^dt 2 + 2e 2S drdt , (26) 

where a,/?, 7,(5 depend on t, r and d£l 2 (9,(f>) is the metric on the 2-dimensional sphere S^ 2 \ Then 
the two-dimensional reduction of the four-dimensional EW theory coupled to a scalar field tp can 
easily be found (here the prime denotes differentiations in r and the dot - in t: 

£( 2 ) = e 2/3 [e- a -''(A 1 - A' ) 2 - e~ a+ ^(^' 2 + i?A\) + e a "^ 2 + ^ 2 A 2 ) - e a+ ^(V + 2A)] + C gr , (27) 

where ip = ip(t,r), V = V(ip), A { = Ai(t,r), A 1 - A' Q = F w and 

C gr = e - Q+2/3 +T(2/3 /2 + 4/3'V) - e a+2 ^(2$ 2 + 40a) + 2ke a+1 (28) 

is the gravitational Lagrangian, up to the omitted total derivatives that do not affect the equations 
of motion. Variations of this Lagrangian give all the equations of motion except one constraint, 

- p' - PP' + dp' + $ 7 ' = ^' + A A l \, (29) 

which should be derived before we omit the 5-term in the metric (taking the limit 5 — > — oo). All 
other equations of motion can be obtained from the effective Lagrangian (|27jf ^l. 

Now, the distinction between static and cosmological solutions is in the dependence of their 
'matter' fields A{ and ip on the space-time coordinates. We call static the solution for which 
Ai = Ai(r) and tp = ip(r). If A = Ai(t) and tp = ip(t) we call the solution cosmological. There 
may exist also the wave-like solutions for which A and tp depend on linear combinations of t and 
r but here we do not discuss this possibility (see, e.g., [25] and references therein). For both static 
and cosmological solutions the gravitational variables in general depend on t and r. To solve the 
equations of motion we may reduce them by separating t and r. It is clear that to separate the 
variables r and t in the metric we should require that 

a = ao(*) + ai(r), p = /%(t) + (3i(r), 7 = 7o (t) + 7l (r), (30) 

Inserting this into the equations of motion one can find the restrictions on the gravitational (and, 
possibly on the matter) variables that must be fulfilled. The details can be found in [21], where 

9 The asymptotic expansions for r — > and r — > oo can be obtained by following the approach proposed below 
for the cosmological solutions. 

10 The equations are equivalent to Eqs.JTTJ) - (|19p but are written in the Einstein frame and in the (t, r) coordinates. 
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one can find the complete list of the static and cosmological spherically symmetric solutions when 
the vector field identically vanishes (we call this case the 'scalar cosmology). Here we only give a 
very brief summary and a simple generalization to nonvanishing vector field. 

The naive cosmological reduction (that supposes all the fields to be independent of r) does not 
give the standard FRW scalar cosmology. As was shown in [20] and [21] (see also the earlier paper 
|15j). all homogeneous isotropic cosmologies should satisfy the following conditions: 

a = $, 7=0, ft" + ke~ 2lSl = , ke~ 2131 - 3ft 2 - 2ft" = C , (31) 

where C is a constant proportional to the 3-curvature (its time dependence is given by the factor 
e~ 2Q °) and the third equation is the isotropy condition. Neglecting inessential constant factors, we 
also have chosen a\ = 71 = 0. We see that for the naive reduction the isotropy conditions in (|3ip 
can be satisfied only if k = and that the first condition is not dictated by (1290 . Thus, the naive 
reduction gives, in general, a homogeneous non- isotropic cosmology. 

For the FRW scalar cosmology 0' ^ and all the conditions (|31|) must be satisfied. Then 
the effective one-dimensional Lagrangian describing both the naive and FRW cosmology is 

£(!) = 6ke a+ ~< - e 2/ V +7 (F + 2A) - e a ~^(2$ 2 + 4$a - ip 2 )] , (32) 

where A; is a new real constant related to k and C; ot,/3, 7 and ip depend only on t. 

Taking now a(t) = (3(t) we get the Lagrangian of the FRW scalar cosmology, for which it is 
not difficult to derive the equations of motion. As in the static case, we can use a as the new 
independent variable and to derive the following first order equation for x(V0 — dip /da = ip'{ot): 



dX 1 (R 2n 



x + 



V + 2A 



(33) 



If we could analytically solve this equation, we would derive the expression for x(V0 an d then, to 
find ip(a), one should solve the equation ip'(a) = x(V'( a ))- The last step must be to write and solve 
the equation for a(t): 

a \t) = [V(iP(a)) + 2A] [6 - ip' 2 (a))- 1 . (34) 

All this looks nice but does not give the desired analytical expressions although the above formulae 
may be useful for a qualitative analysis of some properties of possible solutions. To get approximate 
analytical expressions the approach used for studying the static solutions is much more useful. Of 
course, it can be applied mutatis mutandis to the scalar FRW cosmology but we omit this exercise 
and will write only the solutions for the approximate Lagrangian of the EW model. 

We first write the general cosmological EW Lagrangian supplemented by the minimally coupled 
scalar field (that may represent either matter or inflaton). At first sight, the dimensional reduction 
of the spherically symmetric Lagrangian (j27f )-(|28f) with the vector field must not differ from the 
usual one used for the scalar cosmology and can be written as: 

£« = 6ke a+ ~< + e 2f3 [e- a -^Al - e~ a ^ f A\ - e a+ ^(V + 2A) - e a ^(20 2 + A$a - ip 2 )] , (35) 

where all the fields depend on t. Then, taking a = (3, we apparently obtain a FRW type cosmology 
with the vector field. However, unlike the scalar field, the two-dimensional vecton field equations, 

d [e ao -^+ 2 ^Ax] = _ /1 2 e «0+70+2ft ii ; ( 3 g) 
3^0-70+2/3!^] = _ M 2 e 3a -7o+2/3 Mo _ ( 3? ) 

do give additional constraints on 0i(r). The first equation does not depend on ft, but the 
second one requires either ft(r) = or P[(r) = const. The second condition gives Aq ~ A± and 
so (|29p is incompatible with the isotropy condition a = (3. Therefore, there remains only the first 
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case, P[(r) = 0, from which it follows that k = k = 0. Although the constraint (|29p is identically 
satisfied (as we suppose that 7' = 0) it does not give the necessary isotropy condition a = $ that 
automatically emerges in the scalar cosmology case. As we'll see in a moment, this condition cannot 
be exactly satisfied in the vecton cosmology and can only be approximate^]. 

Summarizing this discussion, we consider the vecton plus scalar cosmology described by the 
Lagrangian (|35|) with k = k = and A\ being the A z component of the four-dimensional vector 
field (it follows that the cosmology must be in general non- isotropic). To write the corresponding 
equations of motion in a most clear and compact form we introduce the temporal notation 

p = h a + 2P), a=\{p-a), A± = e~ 2p+ ^ (A 2 ± p 2 e 2 ~< A 2 ) , V = V(^) + 2K. (38) 
where A\ = A z = A. Then the exact Lagrangian for EW plus scalar cosmology is: 

£ (1) = e 2p-7(^2 _ Qp 2 + g^ 2) + e 3 P -l A _ _ e 3p+ 7 _ (39) 

We see that A, tp, p, a are the dynamical variables and e 7 is the Lagrangian multiplier, variations 
of which give us the remaining energy constraint: 

ip 2 - 6p 2 + 6a 2 + A- + e 27 V = (40) 

(the momentum constraint (|29|) is satisfied by construction). The other equations are: 

i + ( / o + 4o-- 7 )i + e 2 V^ = 0, (41) 

4p + 6p 2 - 4p7 - 6& 2 + \a_ + i? - e 27 V = o , (42) 

a + 3ap-&j- -A- = . (43) 

^ + (3p- 7 )^ + I e 27-^ = 0) (44) 

These equations are much more complex than the equations of the scalar cosmology. They are not 
integrable in any sense and rather difficult for a qualitative analysis. Nevertheless, they are not 
much more difficult than the equations for the static solutions considered in Section 2.2 and we 
may apply the same approach to solving them in asymptotic regions. If the scalar field identically 
vanishes these equations are essentially identical to the equations for the static states. 

They also would be greatly simplified if it were possible to neglect the cr-field. Unfortunately, 
it is evident that this is in general impossible because then A^ = and the last condition is 
incompatible with the other equations. This means that the exact solutions of the EW model 
(even with many scalar fields minimally coupled to gravity) should be non-isotropitJ^ 2 !. In the 
next subsection we consider the simplified model obtained by neglecting the anisotropy. 



2.3.2 Simplified model 

To get the simplified model we neglect Eq. flU]) and take a = 0, ifi = 0, V(ip) = 0. Then p = a = ft 
and the approximate effective Lagrangian (I39|) becomea^l 

C a = -6a 2 e 3a - 7 - 2Ae 3a+7 + iV" 7 - p 2 A 2 e a+ ~< , (45) 

n This fact was not emphasized in the first version of this report. Moreover, the solutions derived there are only 
approximate solutions of the four-dimensional theory as clearly follows from the general equations written below. 
The author is thankful to V. Rubakov for calling his attention to this apparent inconsistency in the presentation of 
the results in the first version. 

12 If one would introduce other scalar fields non-minimally coupled to gravity, this statement may become not 
valid. At this stage of investigation, we are not ready to add other vector fields or fields with the the spin 1/2. 

13 Above we completely neglected the dimensions of all the variables and omitted the gravitational constant. Here 
we only need to restore one of the dimensions supposing that [t~ 2 ] = [k] = [A] = [/i 2 ] = [L~ 2 ]. 
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The corresponding equations of motion are the three equations (|4Q|) - ()42|) with a = ifi = V = and 
p = a. The first equation, (I40p . is equivalent to vanishing of the Hamiltonian. Denoting f = e a 
and taking the gauge fixing condition 7 = (the 'standard' gauge) we havd 1 " 4 ! 

= f [-6 f 2 + 2A f 2 + A 2 + li 2 A 2 } = 0. (46) 

Another useful gauge (the LC gauge) is a = 7. In this gauge, the effective Hamiltonian is: 

H\ = -6/ 2 + 2A/ 4 + A 2 + / 2 p 2 A 2 = . (47) 

It is not difficult to write the equations independent of the gauge choice and we leave this as a 
simple exercise to the reader. 

Let us write the equations of motion in the LC gauge a = 7. In analogy to the static case 
we write them in the first order form (the first equation is the definition of F), 

F = a = f/f, F + F 2 = \kf 2 + \p 2 A 2 , A = B, B = -Li 2 f 2 A, (48) 

6 o 

and the Hamiltonian constraint is a simple polynomial function of /, F,A,B: 

U\ = -6f 2 F 2 + B 2 + 2A/ 4 + Li 2 f 2 A 2 = . (49) 

Similarly to our previous consideration of the static equations, we prefer to change the independent 
variable to a = In/. It is convenient to introduce two new functions, ip(a) and G(a), 

i/j'(a) = F'(a)/F(a), G(a) = F 2 (a) - \he 2a , (50) 
and use the following equations (the prime denotes differentiation w.r.t. a): 

G' + 2G = -li 2 A 2 , A" + ij/A' + fi 2 e 2a F~ 2 A = 0. (51) 
3 

Of course, instead of the first equation we can use the equivalent equation for F 2 = T that directly 
follows from (I48p : 

T' + 2F= \i\e 2a + \p?A 2 . (52) 
Together with the constraint (I49p . rewritten as 

Hi = -e 2a [F(6 - e~ 2a A' 2 ) - 2Ae 2a - [i 2 A 2 \ = , (53) 



equations (|51f) form the complete system describing cosmology in the LC gauge. Note that the 
constraint (|49p is the integral of motion and thus it is sufficient to require that it vanishes just at 
one point, say, at t = or a = —00. To derive possible asymptotic behaviour of the solutions for 
\a\ — > 00 it is natural to expand A in powers of e a and to self consistently use the general solution 
of the first equation, 

G(a) = ^ 2 e~ 2a J daA 2 {a)e 2a , (54) 

with the relations for the expansion coefficients obtained from the second equation. 

In this way we can find, step by step, the asymptotic expansion. In the asymptotic region 
a — ► —00 we can then find the following possible asymptotic behaviour: 

00 ■- 00 -1 00 

A = J2^ne an , T = e~ 2a + ]T f„e m , V' = -1 + £ n^ n e na , (55) 



n=0 



n=2 



n=2 



14 Here we treat cosmological solutions independently of the static states and use somewhat different notation. For 
example, the normalization of A and of fi 2 are slightly different. 
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where Cxm Aq, A\ are arbitrary constants 15 !; A n , T n for n > 2 are derived recursively from (j5T 
(I54p . and ^ n from definition (I50p . The first coefficients are: 



A, 



1 



±\/6, A 2 = 0, A 3 
2 



1 ^ 



V 2 A 



T 2 = -n 2 Al , ^3 = -/j 2 A A 1 ; 



V ; 2 



1 



"^2 ^3 



1 



2 Coo ^Cqq 

Thus we find the differential equation for the metric function /(t) ('scale factor'): 



/C^[f + 2^2/ 2 + 2^ 3 / 3 + -] 



(56) 
(57) 

(58) 



and if we solve it we can find the vector field A(t) by using (|55|) . (|56|) . Neglecting the third term 
in the r.h.s. it is easy to solve this equation finding the dependence of / on t: 



/(t) = ^sinh 



~M 2 ^o ) 



At first sight, the exponential growth of f(t) suggests a possibility of an inflation character of this 
solution. However, this is only the first approximation and we should take into account higher 
order terms to get a more solid conclusion 16 !. Moreover, we see that the qualitative character of the 
solutions essentially depends on the physical parameters Aq, /x 2 on which at the moment we have 
no reliable information 17 !. 

The discussed solution is not unique. Using the above equations we can derive another one, 
for which both A and F are finite for a — > — oo. To get it we take 



G(a) 



daA 2 (a)e 



2a 



(59) 



and then apply the above procedure. Then, using the expansions 



A = ^A n e 

n=0 



2na 



we can find that 



•7"b = F n 



6 



i, 2 A 2 



F=J2 Fne 2na , 

n=0 



Fx = 2F F l 



n=0 



-A + -n 2 A A 1 
6 6 

2 



2na 



A x = -fi 2 A /4F 2 = -3/2A , A 2 = -A^ 2 + 6^i)/16F ( 



: 



(60) 



(61) 



(62) 



where now Aq is the unique arbitrary constant (in the above solution we have one more constant 
Coo). Instead of Eq. (|58p we now have the equation: 



F = a = [Fq + 2F 1 F ' 1 e 2a W + ...]* , 
which can easily be solved in this approximation: 

f = e a = 2e Fot (l - 2F 1 F - 1 e 2Fot y 1 . 



(63) 
(64) 

15 A\ is defined by the constraint (|49[l : putting the first terms of the expansion (|55p into (|49[) or (|53|) we get 
Ai — \/6, which is sufficient for satisfying the constraint. 

16 An interesting exercise couid be to keep four terms in the r.h.s. of (|58|) and express the solution in terms of the 
elliptic functions. The behaviour of f(t) in this approximation essentially depends on all the parameters. 

17 The dependence on A only occurs in the omitted fourth-order terms. 
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The scale factor vanishes if t — > —oo. The parameter F\F 1 strongly depends on Aq, A, and /j, 2 . 
It may be positive or negative, small or large: 

^^(a-^ 2 )^)" 1 - (65) 

We see that for the negative values of FiFq 1 the scale factor / can grow with t up to the maximum 
value l-Fo-F^ 1 ! and then decrease to zero. For positive FiF^ 1 it blows up when the expression 
in the brackets vanishes, but then the approximation is certainly inapplicable. Thus, we must 
be very cautious in making definite conclusions basing on this simple result. The approximation 
(I63p can only be reasonable when |/ 2 < Fq. As Fq and F\ strongly depend on A, on the 
absolutely unknown mass [i and on the arbitrary constant Aq, it is not possible (at the moment) 
to make conclusive statements on the general properties of this solution though it depends on less 
parameters than the first one. Note only that both solutions are compatible with existence of a 
period of fast growing scale factor. 

2.3.3 A generalization 

The expansions for both cosmological solutions can generally be written as follows. Let us use the 
more convenient notation: 

oo n 

F{a) = F 2 (a) = C^eT 2a + ]T T n eT na , A n = ^A x A n ^ , A = A/3. (66) 

n=0 1=0 

Then the general expression for T n can be written as 

T n = ii 2 An/ [3(n + 2)] - kd n>0 + AS n , 2 , (67) 
and Eq. (l5ip for A(a) = ^2, n A n e aa rewritten as 

2FA" + T'A 1 + 2fi 2 e 2a A = , (68) 
defines the recurrence relation for n > (A n = for n = —1, —2): 

n 

2Coo (n + 1) (n + 2) A n+2 + ^ m (m + n) T n - m A m + 2 / u 2 A n _ 2 = 0. (69) 

m=l 

Taking n = we find that Cqo^ = 0. The first solution discussed above is obtained when A 2 = 0, 
the second - when = (in the second case T 2 n+\ = A 2n +i = for all n > 0). In both cases Aq 
is arbitrary and other A n are derived recursively from Eqs. ([67j) . (f69l) . taking into account that for 
the first solution A\ is defined by the constraint (f5"3"|) . i.e. A\ = i-^/6. 

Note that Coo > but the signs of Aq and of A\ are not fixed. For this reason, the structure of 
the series expansion of A(a) is rather complex. It follows that the behavior of a(t) with growing 
time may change its character (for example, depending on the parameters Coo, Aq, fi 2 , A, it may 
change the exponential growth to oscillating or even chaotic behavior). At the moment, the main 
unsolved problem is to derive the asymptotic behavior of A n and T n for n — > oo. An educated 
guess is that the main terms are given by some powers of n and thus the expansions in powers of 
f = e a have a finite radius of convergence in the /-plan^f). If this is true, the radius of convergence 
and the precise position and nature of the corresponding singularity in the complex /-plane could 
give us a very important information on the solutions. 

18 This was proved for the analogous expansions near the horizons, see [22] 
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The dimensionless functions J- (a) / p? and A(a) of general solution depend on the dimensionless 
constants A = A/// 2 , p = C^/ p 2 and Aq. The dependence on Aq is especially simple for the second 
solution. It is not difficult to prove (e.g., by induction) that 



A 2 k/A = a 2 kA 2k , T )Tq = f 2 kA 2k , Fq 



where a 2 k, jik depend only on A. We thus have the expansions 

oo 

A(a)/A = l + Y / a2ke 2k «, TjT Q 

k=l 



l + E hk < 

k=l 



2ka 



a 



a 



In Ar 



in which a 2 k, hk can De derived by the recurrence relations 



(70) 



(71) 



hn = o-^l" + E( n + X ) la 2ka 2n -2k, 
6 k=0 



E K n + k)f 2n ~2k a 2k + 3a 2n -2 = . 



k=l 



It follows that the equation for a(t) can be rewritten as 



da 
~d7 



1 + E hk< 

k=l 



2ka(r) 



T = p Aq I y/Q t . 



(72) 



(73) 



An interesting property of these expansions is the following: if the dimensionless parameter Aq is 
large enough, e a may remain small even for large positive a (thus the scale parameter e a may be 
very large). This means that it may be possible to approximate the expansions by a small number 
of terms (in some domain of t in which e a is very large). This property is of interest for discussing 
inflation models. The first several terms in the expansions are easy to derive (A = A//^ 2 ): 



a 2 



3 

2 ' 
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d4 



A 



h 



f - 3 
/4 "4 



a 6 



-A 



A 



21 
20 



A -20 



(74) 



(75) 



Our simplified Lagrangian (|45p should be considered as a (mainly) mathematical model for a 
preliminary study of the EW theory. To really discuss its cosmological applications we must first 
find an approximation valid for high values of / and then return to the complete set of the equations 
of motion, to which the simplified model is only a rather crude approximation. We also should not 
forget that it is absolutely necessary to include into consideration 'ordinary' matter before one can 
really discuss physical picture of the cosmological evolution. 



3 Discussion 

In this paper we briefly summarized the main ideas of the Einstein - Weyl model and presented 
its new interpretation, as well as some results obtained investigating its simplest solutions. We 
only considered the static spherically symmetric solutions and, in cosmology, only an artificially 
simplified homogeneous model. As we noted in [21], even small deviations from the spherical sym- 
metry may result in a qualitatively different theory. In particular, if we consider axially symmetric 
configurations infinitesimally deviating from the spherically symmetric ones, we will find additional 
scalar fields in the vecton dilaton gravity, which may be very important in cosmological considera- 
tions and in analysing black holes. We did not touch these problems here. Moreover, even in the 
spherically symmetric case our study is incomplete. In the static case, we have only proven that 
there may exist two horizons and derived the solutions close to the horizons. In cosmology, we have 
studied only the asymptotic behavior of the solutions in the simplified model. 
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As we mentioned above, we expect that the complete solutions should reveal some sort of chaotic 
behavior. To study these phenomena we must first carefully discuss the physical parameters of the 
theory. In the original formulation these are: the gravitational constant, the cosmological constant, 
and the vecton mass. In addition, the asymptotic boundary conditions introduce other parameters, 
the dependence on which is highly nontrivial. This does not allow us to make sound conclusions (or, 
even guesses) about the global behavior of the solutions derived in our essentially local approach. 
For example, if we try to glue together the left and the right asymptotic approximations, we will find 
that the gluing procedure is strongly dependent on the parameters that characterize the influence 
of the nonlinear terms in the equations, up to producing chaotic effects. This requires a very careful 
qualitative and numerical study of the equations. Of course, the most important task is taking into 
account the 'ordinary' matter. 

Finally, we must admit that the vecton field is a rather unusual feature of the Einstein- Weyl 
model. I have found just a few papers in which a massive vector field is introduced in the cosmo- 
logical context (see, e.g., |26|-|28| and references therein"*^!) . Thus the unified model of dark energy 
and vecton dark matter considered here looks as fresh and new as it was in 1923. In addition, I wish 
to stress that the original Einstein Lagrangian (j3]) or (JHJ) is more interesting and exciting than the 
simplified theory (|12|) (in particular, one may expect that in the original formulation there exists a 
relation among dimensional parameters that are arbitrary in the theory (1120 ). Unfortunately, the 
original theory is much more difficult to deal with and thus the prime goal must be the study of 
the simplified theory. In this paper I only give a sketch of how to begin such a study. 

Additional comments 

1. To better understand the brief exposure of the static solutions in Section 2.2 the reader is 
advised to consult our paper [22J. There we proved that the expansions of the solutions represent 
analytic functions that analytically depend on all the parameters. This means that the two horizons 
defined by Eq. ([25j) and corresponding to the same parameters Aq, Bq, Fq, belong to the same static 
solution in spite of the fact that it is represented by two different expansions near the two horizons. 
Also, analyticity tells us that any solution can be analytically continued from one horizon to the 
other as well as to any regular point in the interval < r < oo. The same remarks are also relevant 
to the cosmological solutions considered in Section 2.3. Recalling that by crossing the horizon 
we pass from the static to cosmological solutions one may even use the static coordinates for and 
alternative description of cosmologies (see, e.g., [29], [30]). 

2. Recently, several authors considered massive vector models of inflation (see, e.g., |31j-|34j 
and references therein). As distinct from the EW model and from the first paper [27] on the vector 
inflation, they introduce some non- minimal couplings of the vecton to gravity and/or many vector 
fields. In the Einstein approach additional vector fields (as well as the inflaton - type scalar fields) 
can possibly be produced in some higher-dimensional version. It is not clear whether non-minimal 
couplings can be produced but, possibly, similar effects could be imitated in the really non-linear 
Eddington-Einstein type theory. There exist attempts to use the non-linear Eddington-Einstein 
type actions for description of dark energy and dark matter (see, e.g., [SS] and references therein). 
Note also that, possibly, the anisotropic effects of the vector field can be smeared out by inflation 
(see, e.g., [SB] where the evolution of the universe with the primordial magnetic field was discussed). 

3. In this paper, we did not discuss the possible interpretation of the vecton as a dark matter 
candidate. One reason for this is our inability to estimate its mass, even by order of magnitude. 
The other reason is that the vecton is interacting only with gravity and thus can be produced (in 
abundance) only in very high gravitational fields. This requires essentially quantum considerations 
that are beyond the scope of the present paper. 

19 In the second version of this paper I omitted two references in which the massive vector has only Ao(t) component 
and thus Fm = 0. 
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4 Appendix 



The connection ([5]) is a special case of the general expression for the connection in affine spaces. 
The most general symmetric affine connection has the form (see, e.g., |37j ) : 

T kl = \[s mn {Snk,l + Sln,k ~ S k l,n) + S mn (s nM + S lnk - S Hn )] , (76) 

where Ski is an arbitrary symmetric tensor, s mn is the inverse matrix to sm, and Skin is an arbitrary 
tensor that is symmetric in k and I. Both the Weyl and Einstein connections belong to the subclass 
for which Skin can be presented in the form: 

Skin = as k ii n + PiSnkU + Slnik) ■ (77) 

We may call it the Weyl-Einstein connection (defining the Weyl-Einstein spaces). 
Inserting (f77|) into (fT6"|) . we find: 

T Ti = \ K n (snk,i + s lnik - s kl>n ) + a(C il + SP i fc ) - (a - 2(3)s kl i m ] . (78) 

Now it is easy to find that the Einstein connection ([5]) corresponds to a = —(3 = ^. The Weyl 
connection introduced in [S] corresponds to a = 1, (3 = 0. I could not find a discussion of the 
geometry of spaces with the Einstein connection in accessible literature. The geometry of the 
Weyl spaces is considered in [S] {D = 4, Lorentzian signature), and in [37] (D = 2,3,4, Euclidean 
signature). 
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